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The proton ordering model for the KH 2 PO 4 type ferroelectrics is modified by taking into account the depen¬ 
dence of the effective dipole moments on the proton ordering parameter. Within the four-particle cluster ap¬ 
proximation we calculate the crystal polarization and explore the electrocaloric effect. Smearing of the fer¬ 
roelectric phase transition by a longitudinal electric field is described. A good agreement with experiment is 
obtained. 
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1. Introduction 

The electrocaloric (EC) effect is the change of temperature of a dielectric at an adiabatic change of the 
applied electric field. Research in this field is driven by a quest for materials that can be used for efficient, 
environment-friendly, and compact (on-chip) solid-state cooling devices. 

The current state of the art on the electrocaloric effect research for ferroelectrics is well summarized 
in fll. At the moment, the largest effect is observed in perovskite ferroelectrics. Thus, in in the 
PbZr 0 . 95 Ti 0 . 05 O 3 thin film with a thickness of 350 nm in a strong electric field (480 kV/cm) the obtained 
electrocaloric temperature change is AT = 12 K. Ab initio molecular dynamics calculations (3l predict 
AT » 20 K in LiNbOs. In the hydrogen bonded ferroelectrics of the KH 2 PO 4 (KDP) type, the electrocaloric 
effect was studied for relatively low fields only. Thus, it has been obtained that AT 0.04 K at i?» 4 kV/cm 
dl, AT a: 1 K at £ r: 12 kV/cm 0, and AT r: 0.25 K at Tc and E » 1.2 kV/cm 0]. 

Theoretical calculations of the electrocaloric effect in KDP have been made in 0 within the Slater 
model d and In the paraelectric phase only. It is also known that the Slater model gives incorrect results 
in the ferroelectric phase, and more complicated versions of the proton ordering model are required for 
an adequate description of these crystals. Thus, the effect of electric field on the physical characteristics 
of the KDP type crystals, such as polarization, dielectric permittivity, piezoelectric coefficients, elastic con¬ 
stants, has been described within the proton ordering model with the piezoelectric coupling to the shear 
strain ee floUlzIl and with proton tunneling Q] taken into account. However, these theories required, 
in particular, invoking^two different values of the effective dipole moments for the paraelectric and fer¬ 
roelectric phase (^[ 1 ^. This made impossible a correct description of the system behavior in the fields 
high enough to smear out the first order phase transition. There is an inner logical contradiction in the 
model: whUe no physical characteristic of a crystal should exhibit any discontinuity in the fields above 
the critical one, there is no smooth transition between the values of model parameters, rigidly set to be 
different for the two phases. 

In the present paper we suggest a way to remove this contradiction. Assuming that the difference 
between the dipole moments is caused by non-zero values of the order parameter, we modify the proton 
ordering model accordingly. The field dependences of polarization, smearing of the first order phase 
transition, and the electrocaloric effect are described. 
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2 . Thermodynamic characteristics 

We consider the KDP type ferroelectrics in the presence of an external electric field £3 applied along 
the crystallographic axis c, inducing the strain ee and polarization P 3 . The total model Hamiltonian reads 

H^NHo + Hs, ( 1 ) 


where N is the total number of primitive cells. The “seed” energy Hq corresponds to the sublattice of 
heavy ions and does not exphcitly depend on the proton subsystem configuration. It is expressed in terms 
of the strain ee and electric field £3 and includes the elastic, piezoelectric, and dielectric contributions (ijl 

Ho = - e36^3£6 - , ( 2 ) 

where v is the primitive ceU volume; e^g, are the “seed” elastic constant, piezoelectric coefficient, 
and dielectric susceptibility, respectively. 

The pseudospin part of the Hamiltonian reads 

E ///'W^7')^^ + 4h + E2V'6£6^-EM/-B3^ + H£. (3) 

Z rt Z Z f Z f Z 


Here, the first term describes the effective long-range interactions between protons, including also in¬ 
direct lattice-mediated interactions (3 Ell; cFqf is the operator of the z-component of a pseudospin, 
corresponding to the proton on the /-th hydrogen bond (/ = 1, 2, 3, 4) in the ^-th cell. Its eigenvalues 
Gqf = ±1 are assigned to two equilibrium positions of a proton on this bond. 

In (1^, Hsh is the Hamiltonian of short-range interactions between protons, which includes terms lin¬ 
ear over the strain jn]] 


Hsh = 


E { ( {(^ql + Crq2 + CTqO + qi) 


(Ss 5i 1 

+ I — ^6 - "^^6 I {(Tql(Tq2(^qO + ql(^q2(Tqi + (^ql(Tq3(^qi + (^q2(^qO(TqA) 

+ “(i^ + 8a£o){CTq\(^q2 + CTq^O'qA) + -{V - 8a£6)i0'q20'q3 + O'qA^'ql) 


U (p 1 

+ — iOqi(Tq3 + aq2(Tq4) + —(Tq\aq2(Jq3aq4j. 


(4) 


Here, 

1 1 

U-—WI—E, ^ - 4e- 8w+ 2wi, 

and E, w, wi are the energies of proton configurations. 

The third term in (3) is a linear over the shear strain ee field due to the piezoelectric coupling; y/e is the 
deformational potential. The fourth term effectively describes the system interaction with the external 
electric field £ 3 . Here, [if is the effective dipole moment of the /-the hydrogen bond, and 


£1 = A «2 = A«3 = £4 = £■ 


The fifth term in is introduced in the present paper for the first time. It takes into account the 
assumed dependence of the effective dipole moment on the order parameter (pseudospin mean value) 


He = 



2 



2 


(S) 


It is equivalent to a term proportional to £|£3 in a phenomenological thermodynamic potential. Note that 
the terms hke £f £3 are not allowed because of the symmetry considerations, and we keep the Hamilto¬ 
nian to be linear in the field £ 3 . 
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In view of the crystal structure of the KD P type f erroelectrics, the four-particle cluster approximation 
is most suitable for short-range interactions lisilldl . Long-range interactions and the term He are taken 
into account in the mean field approximation. Thus, 


He = -UN^i'Estj' 


/=! 


16N^i'E3T]^ 


(6) 


Combining the fourth term in and the first term in ( 6 ), we obtain the following term in the Hamiltonian 
-(/i-h l 2 fi'rf']E 3 Y.qf(Tqfl 2 . Effectively, the term 12[j.'t]^ in {fi+ \2^'Tf) describes the jump of the dipole 
moment at the first order phase transition, its different values for the paraelectric and ferroelectric phase, 
and its smooth behavior in the fields above the critical one, when there is no jump of 77 . We can now use 
a single value of /i for both phases and remove the logical contradiction of the earlier theories, described 
in Introduction. 

Proceeding with the standard calculations of the cluster approximation (l^ EE EEl. we obtain the 
following expression for the proton ordering parameter 


r] = {(Tql) = {(Tqz) = (CTqs) = (crqi) = — , 


where 


m - sinh( 2 z-i- pdsEe) -i- 2 hsinh(z-)Sdie 6 ). 

D - cosh(2z + pdsEe) + 4hcoshfz - fid 1 ee) + 2flcosh pSaEe + d, 

1 1-1-77 Bu. , n 

z = - In -—- + pvcT] - Py/eee -t ^£3 -t 676/7 t]^E 3 , 

2 1 - 77 2 


4vc = /ii(0) - 1 - 2 / 12 ( 0 ) -H /i 3 ( 0 ) is the eigenvalue of the long-range interactions matrix Fourier transform 
///'=lR,-R,,///'(^^7');76 = i/fcBr. 

The thermodynamic potential is then obtained in the following form 

G = |c^6°£6-«'«^36e6£3-|li3^|+2Vc772-Hl6/7%77^ (7) 

2 2 2 
-H-ln2--In(1 - 77^)--InD - ncTeee ■ 

Here, ae is the formally introduced shear stress conjugate to the strain eg. In numerical calculations we 
put (jQ - 0. The condition of the thermodynamic potential minimum 


yields an equation for the strain eq 


f-1 

W^e/r.Ea.o-e 


Fo D 

CT6 = C66 £6-e^eEs+ —!] + — . 


(8) 


In the same way, we derive the expressions for polarization P 3 and molar entropy of the proton subsys¬ 
tem 


P 3 


S 


--f-1 

V \dE3)T,ae 



2 \dTjE,,a, 


2 M 

-ln 2 -i-ln(l - T]^) +1x10 + 2 TzTr] -H — 


Here, A/a is the Avogadro number; R is the gas constant. The following notations are used: 


(9) 

( 10 ) 


r^-dsMs-6aMa+diMi, 

1 , 2 
zr = --p^(Vc77-7/76e6 + 6/7 77 £ 3 ), 
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M- Abpwcosh[z- pSIee) + pwid + 2apecoshpda£6 + P^er, 

Ma - 2ash\hfi8a£B,Ms - smh[2z + fi8s£s),Mi - 4i)sinh(z - 

Expressions for dielectric susceptibilities, piezoelectric coefficients, and elastic constants derived 
from equations (8), (9) are slightly different from the previous ones (3, where the dependence of the 
effective dipole moment on the order parameter was not taken into account. Numerical calculations, 
however, showed (3 that in zero electric field fhe difference is minor. 

The molar specific heat of the subsystem described by the Hamiltonian JT) is 

AC" = r(^) =r(Sr + S^pr + S,£r). (11) 

Here, 


f5S) R 



2TzTiq6-qM) + Ne- — 


= — [DTzT + {qs-'nM)Zr,], 

Notations introduced here are described in appendix. 

Then, the total specific heat is 

C^AC^ + Cregular- (13) 

Here, is assumed to describe all the anomalies of the specific heat at the phase transition, whereas 
the regular background contribution to the specific heaf, mostly from the lattice of heavy ions, is approx¬ 
imated by a linear temperature dependence 


, , M 

-2[qe-r]M)yfe-:i+ —r 


( 12 ) 



Cregular — Co + Cl (T' 7c)- 


(14) 


As wih be discussed later, this hnear approximation agrees with the experimental data. 
Finahy, the electrocaloric temperature change is calculated using the known formula 



(IS) 


where the pyroelectric coefficient is 


dPs] 0 2{iu+12i^'q^) 


(16) 


V - vNp,!2 is the molar volume. 


3. Numerical calculations 

To perform the numerical calculations we need to set the values of the following theory parameters: 

— the Slater energies £, w, wi; 

— the parameter of the long-range interactions vd 

— the effective dipole moment p and the correction is due to proton ordering /i'; 

— the deformation potentials iffs, 8s, 8a, 5i; 
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— the “seed” dielectric susceptibility X 33 , elastic constant Cgg, piezoelectric coefficient e^gj 

— the parameters of the lattice specific heat Co and Ci. 

They are chosen, obviously, by fitting the theoretical thermodynamic characteristics to the experimental 
data, as described in (3- The obtained optimum sets of the model parameters are given in table[l] 

To describe crystals with different deuteration levels, we use the mean crystal approximation, where 
the theory parameters are assumed to be linearly dependent on deuteron concentration (except for the 
parameter Vc, for which a small deviation from the linear dependence is assumed, as it is chosen from 
the condition that the calculated transition temperature coincides with the experimental one, which is 
also slightly non-linear). The dependence of the energy levels and interparticle interaction constants on 
deuteration is caused by the corresponding geometrical changes in the crystal structure with deuteration 
(elongation of the hydrogen bonds, changes in the distance between the equilibrium positions of H or D 
on the bonds, changes in the lattice constants, etc). 


Table 1. The optimum sets of the model parameters for different crystals. As KD 2 PO 4 we denoted 
K(Hi_j;D;c) 2 P 04 with X = 0.89. 



npU 

(K) 

£/A;b 

(K) 

w/ks 

(K) 

Vc/fcfi 

(K) 

(10-3“ 

m' 

(10-3° c-m) 

I33 

KH 2 P 04 

122.22 

56.00 

430.0 

17.55 

5.6 

-0.217 

0.75 

KD 2 P 04 

211.73 

85.33 

730.4 

39.26 

6.8 

-0.217 

0.39 

KH 2 AS 04 

97 

35.50 

385.0 

17.43 

5.5 

-0.033 

0.7 

KD 2 AS 04 

162 

56.00 

690.0 

31.72 

7.3 

- 0.000 

0.5 



y/elks 

(K) 

Sslks 

(K) 

Salk^ 

(K) 

S\lk^ 

(K) 

c.h'D 

‘'66 

(10° N/m^) 

e° 

‘^36 

(C/m3) 

Co 

J/(mol K) 

Cl 

J/(mol k3) 

KH 2 PO 4 

-150.00 

82.00 

-500.00 

-400.0 

7.00 

0.0033 

60 

0.32 

KD 2 PO 4 

-139.89 

48.64 

-1005.68 

-400.0 

6.39 

0.0033 

93 

0.32 

KH 2 ASO 4 

-170.00 

130.00 

-500.0 

-500.0 

7.50 

0.01 

60 

0.32 

KD 2 ASO 4 

-160.00 

120.00 

-800.0 

-500.0 

6.95 

0.01 

98 

0.40 


The primitive ceU volume is taken to be n = 0.1946 ■ 10“^^ cm^ for K(Hi_jcDj:) 2 P 04 and v - 0.202 ■ 
10“^^ cm^ for K(Hi_j:Dj:) 2 As 04 , irrespectively of the deuteration. The energy wi of proton configurations 
with four or zero protons near the given oxygen tetrahedron should be much higher than e and w. 
Therefore, we take wi -00 [d-0). 

As we have already mentioned, when the dependence of the effective dipole moment on the order 
parameter is taken into account, the agreement between the theory and experiment for most of the cal¬ 
culated dielectric, piezoelectric, elastic characteristics, and specific heat of the studied crystals in the 
absence of an external electric field is neither improved nor worsened (see El). However, the present 
model allows us to describe more consistently the smearing of the first order phase in high electric fields. 

The temperature dependence of the specific heat of KH 2 PO 4 and KD 2 PO 4 is shown in figure[l] The con¬ 
tribution AC" is essential in the transition region and satisfactorily describes the experimental anoma¬ 
lies. As one can see, the total specific heat above Tc can be well approximated by a linear temperature 
dependence, thus justifying the linear dependence of Creguiar. given by equation JT4l . 

In figures|2]and[3]we plotted the temperature variation of polarization of K(Hi_;cD;c) 2 P 04 in different 
fields. The agreement with experiment is better at x = 0.89 (and 0.84, see El) than at x = 0 . We believe 
this is due to proton tunnelling, essential in non-deuterated samples, which is not included in our model. 

The field E 3 , which in these crystals is the field conjugate to the order parameter, induces non-zero 
polarization P 3 above the transition point. Polarization has a jump at Tc, indicating the first order phase 
transition. With an increasing field, the polarization jump decreases, whereas the transition temperature 
Tc increases almost linearly. The corresponding dTcldE^ slopes are 0.192 and O.llS K cm/kV for x = 0 
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Figure 1. The temperature dependence of the molar specific heat of K(Hi_;tDj;) 2 P 04 at x = 0.0 — o (H, 

□ (19)]; at X = 0.86 — A (^. Dashed lines 1’ and 2’: the theoretical results of (inl¬ 
and X - 0.89, respectively (c.f. 0.22 and 0.13 K cm/kV from our earlier calculations and experimen¬ 
tal 0.125 K cm/kV of for x = 0.89). At some critical field E*, the jump vanishes, and the transition 
smears out. The calculated coordinates of the critical point are E* - 125 V/cm, =122.244 K for x = 0 
and 7.1 kV/cm, 212.33 K for x = 0.89, which agrees well with the experiment K should be noted 

that in our previous calculations Gi it was impossible to obtain a correct description of the polariza¬ 
tion behavior in the fields above the critical one, due to the necessity of using two different values of the 
effective dipole moment p in calculations. 

The calculated electrocaloric changes of temperature AT of the K(Hi_;cD;c) 2 P 04 and K(Hi_xDj;) 2 As 04 
crystals with the adiabaticaUy applied electric field are shown in figures|4]and[S] The experimental data 
of 01 were obtained at T = 121 K, which was very close to the transition temperature of the sample used 
in the measurements. 



Figure 2. The temperature dependence of polariza¬ 
tion of KH 2 PO 4 at different i: 3 (MV/m): 0.0 — 1, a 
d; 0.581 — 2, O (^; 1.250 — 3, □ 2.031 — 4, 

0 ( 2 ^. Symbols are experimental points; solid lines: 
the present theory; dashed Unes: the theoretical re¬ 
sults of Ilih. 



Figure 3. The temperature dependence of polariza¬ 
tion of K(Hi_;cDj:) 2 P 04 at x = 0.89 and at different 
E 3 (MV/m): 0.0 — 1; 0.282 — 2, o; 0.564 — 3, □; 0.71 — 
4; 0.846 — 5, 0; 1.128 — 6, A. Symbols are experimen¬ 
tal points taken from (i^l; lines: the present theory. 
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Figure 4. The field dependence of the electrocaloric temperature change of K(Hi_jcDjc) 2 P 04 for x = 0.0 
(soUd Unes) and jc = 0.89 (dashed lines) at T - T° = -2.04 K — 1, l', □; T = 7° — 2, 2', o; r - r° = 3.2 K — 
3, 3 ', o. Experimental points are taken from @| — o, □ and (3l — o. 


As one can see, at small fields (figures|4][^ left-hand) the calculated electrocaloric temperature change 
is a linear function of the field below (curves 1, 1') and a quadratic function above (at 2, 2'). The 
experimental behavior below is not linear at fjs < 2 kV/cm due to the domains: The domains, whose 
polarization is oriented along the field, are heated, whereas the domains, polarized in the opposite di¬ 
rection are cooled, thus the resulting net change of the sample temperature is close to zero. The experi¬ 
mental data for the electrocaloric temperature change at and above T® available for KH 2 PO 4 , as well as 
the ATlAE ratio below at fields above 2 kV/cm (when the sample is in a single-domain state), are well 
reproduced by the theory. 

At higher fields (figures S) [S] right-hand) the calculated electrocaloric temperature changes at tem¬ 
peratures above are larger than below T^. The obtained curves deviate from linear and quadratic 
behavior and reach saturation at E » 500 kV/cm. It should be mentioned, however, that these curves are 
calculated with the linear over the field E 3 pseudospin Hamiltonian J3). It would be very interesting to 
compare our results at high fields with experiment, for instance, to find out when non-linear contribu¬ 
tions to the Hamiltonian cannot be omitted any longer. Unfortunately, no experimental data for AT in 




Figure 5. The field dependence of the electrocaloric temperature change of KH 2 ASO 4 (soUd Unes) and 
KD 2 ASO 4 (dashed Unes) at T-T° = -2.04 K—1, l'; r= T° —2,2'; T- T° =3.2K — 3, 3'. 
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Figure 6. The temperature dependence of the electrocaloric temperature change of K(Hi_xDx) 2 P 04 for 
X = 0.0 (left-hand) and x = 0.89 (right-hand) in different fields. 


the fields above 1 kV/cm are available. And, of course, possibilities for experimental measurements are 
limited by the dielectric strength of the samples. 

As one can see from the temperature dependence of AT (figure [6) for K(Hi_j;Dx) 2 P 04 crystals, the 
calculated electrocaloric temperature change is the largest at temperatures below but close Tc and 
can exceed 6 K; however, the fields required to reach AT that high are not accessible in reality, because 
most likely they exceed the dielectric strength of the crystals. 


4. Conclusions 

Taking into account the dependence of the effective dipole moment on the order parameter within 
the framework of the proton ordering model allows us to correctly describe the smearing of the ferroelec¬ 
tric phase transition in high electric fields as well as the electrocaloric effect in the KDP family crystals. 
The theory predicts the values of the electrocaloric temperature change of a few Kelvins in high fields. 
Additional experimental measurements of AT in the fields above 2 kV/cm are necessary. 


Appendix 

The notations introduced in equations are as follows: 

Ns - 2a{pe)^coshpSa£6 + ^b[pw)^cosh[z- pSiEs) + {Pwi)^d + 2£s(f'(-£8aMa + iv8\M\) 
+ Eg [2a((i8af' coshp8a£6 + [fi8sf cosh(2z -H fi8s£s) + 4h(/l(5i)^ cosh(z - jS^iEg)], 

qs = 2bpwsmh{z- pSi£s)+£sl3 [-<5s coshi2z + (38gEs) + 2bSi cosh(z-;65ie6)], 

A - -pEdaMa + Plu8\M\ +EsP COSh(2z-l- P8sEs) +2a8^ coshpSaE6 + ^b8\ COSh(z-/i5i£6)] , 


riT^Pe + 


2{^i+\2^i'q^) 


(636 - eqgler. 


et- 


In turn. 


vD\ ■’ D V V D 


£ I 2xTzt+ lqs-qAd] 
“ T D-2xz„ 


‘^ 66 - 
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c£. is the isothermal elastic constant at a constant field 




^ 8y/e + fe) 


V D-2zjjX vD[D -2zriX) 



and 636 is the isothermal piezoelectric coefficient 


636 = -[^l 

ldi?3jr,ce \deejT,Es 


0 2(At+12AtV) pee 

V D- 2Zri>{ 


with 


06 = -2x\f/e + fe, fe - 6sCosh{2z + pdsEe) -2b6\ cosh(z-/i5i£6) 
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EiieKTpoKayiopi/mHi/iM e(|>eKT y Kpi/icrayiax Ti/iny KH2PO4 

A.C. BflOBi/i'-P, A.n. MoTheP, P.P. JleBi/mbKniP, I.P. SaneiP 

^ iHCTi/iTyr c|)i3HKM KOHfleHCOBaHMX ci/icreM HAH yKpaiHH, By;i. I. CBEHpipbKoro, 1, 79011 /IbbIb, VKpaiHa 
^ HapiOHa/ibHUM yHiBepcuTer "/IbBiBCbKa no/iirexHiKa", By;i. C. BaHflepw, 12, 79013 /IbBiB, YKpaiHa 

B MOfleyii npoTOHHOro BnopaflKyBaHHfl fl/ia Kpucra/iiB Twny KH2PO4 BpaxoBano 3 a;ie>KHiCTb ec|)eKTHBHHx flw- 
no/ibHMx MOMeHTiB Bifl napaMCTpa npoTOHHoro BnopflflKyBaHHa. B Ha6;in>KeHHi HOTupuHacTHHKOBoro K/iacre- 
pa poapaxoBaHO no/iapwaapiK) Kpucra/iiB ra floc/iiflxeHO e/ieKrpoKayiopuHHMM eefieKT y hux. OnucaHO po 3 MM- 
BaHHfl cerHeroe/ieKTpMHHoro c|)a 30 Boro nepexofly noaflOBXHUM e/ieKipuHHUM no;ieM. OrpuMaHO flo6pe yaro- 
flxeHHfl 3 eKcnepi/iMeHTa/ibHMMu flanuMn. 

KjiioMOBi c;ioBa: ejieKTpoKanopi^>^Htnii ecjieKT, KDP, Knacrepne Ha6yimeHHa, nompi/iaapm 
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